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ABSTRACT 

We present a simple analytic model for the various contributions to the non-thermal emission from shell 
type SNRs, and show that this model's results reproduce well the results of previous detailed calculations. 
We show that the > 1 TeV gamma ray emission from the shell type SNRs RX J1713.7-3946 and RX J0852.0- 
4622 is dominated by inverse-Compton scattering of CMB photons (and possibly infra-red ambient photons) by 
accelerated electrons. Pion decay (due to proton-proton collisions) is shown to account for only a small fraction, 
< 10~ 2 , of the observed flux, as assuming a larger fractional contribution would imply nonthermal radio and 
X-ray synchrotron emission and thermal X-ray Bremsstrahlung emission that far exceed the observed radio and 
X-ray fluxes. Models where pion decay dominates the > 1 TeV flux avoid the implied excessive synchrotron 
emission (but not the implied excessive thermal X-ray Bremsstrahlung emission) by assuming an extremely 
low efficiency of electron acceleration, K ep < 10 -4 (K ep is the ratio of the number of accelerated electrons and 
the number of accelerated protons at a given energy). We argue that observations of SNRs in nearby galaxies 
imply a lower limit of K ep > 10~ 3 , and thus rule out K ep values < 10~ 4 (assuming that SNRs share a common 
typical value of K ep ). It is suggested that SNRs with strong thermal X-ray emission, rather than strong non- 
thermal X-ray emission, are more suitable candidates for searches of gamma rays and neutrinos resulting from 
proton-proton collisions. In particular, it is shown that the neutrino flux from the SNRs above is probably too 
low to be detected by current and planned neutrino observatories. Finally, we note that the magnetic field value 
implied by the comparison of X-ray to gamma-ray emission, ~ 10 /iG, can be used to constrain magnetic field 
amplification. 

Subject headings: gamma-rays: theory — supernova remnants: general: individual RX J 17 13.7-3946, RX 
J0852. 0-4622 — cosmic rays — neutrinos 



1. INTRODUCTION 

For a long time it is believed that the galactic cosmic rays 
observed up to the 'knee' energy (~ 10 ir eV) are accelerated 
in supernova remnants (SNRs, see e.g. Axford 1994). Therel- 
ativistic protons (and electrons) are believed to be accelerated 
by the diffusive (Fermi) shock acceleration (DSA) mecha- 
nism (for reviews see Drury 1983; Blandford & Eichler 1987; 
Malkov & Drury 2001). Strong evidence for electron accel- 
eration to high energies in SNRs was established by observa- 
tions of non-thermal X-ray emission which was attributed to 
synchrotron radiation of multi TeV electrons (Koyama et al. 
1995), while direct evidence for ion acceleration in SNRs has 
not been presented so far. 

Recently, unambiguous detection of > 1 TeV 7-rays has 
been made from the shell-type SNRs RX J1713.7-3946 (Mu- 
raishi et al. 2000; Enomoto et al. 2002; Aharonian et al. 2004, 
2006) and RX J0852.0-4622 (Katagiri et al. 2005; Aharonian 
et al. 2005, 2007), providing the first direct proof for the ac- 
celeration of particles at SNRs to multi- TeV energies. There 
are two candidate emission processes that can account for this 
radiation, namely inverse Compton (IC) of radio and infra red 
photons by multi TeV accelerated electrons or pion decay as a 
consequence of proton-proton (PP) interactions of multi TeV 
accelerated protons with ambient target protons (e.g. Drury et 
al. 1994). 

The distinction between the two mechanisms has impor- 
tant consequences for understanding particle acceleration and 
magnetic field amplification in SNRs. If it would turn out that 
the source is PP emission, this would be the first direct evi- 
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dence for proton acceleration in SNRs. An IC source would 
allow a rather accurate estimate of the downstream magnetic 
field value by comparing the X-ray to gamma ray fluxes (see 
e.g. Aharonian et al. 2006). 

Broad-band emission models with different levels of so- 
phistication where applied in order to analyze the observed 
non-thermal radiation from these SNRs, reaching different 
conclusions as to the dominant > 1 TeV 7-ray emission 
mechanism. For RX J1713.7-3946 Aharonian et al. (2006), 
Berezhko & Volk (2006) and Moraitis & Mastichiadis (2007) 
claimed that PP emission is favorable and IC is unlikely, 
Porter et al. (2006) claimed that IC emission is consistent. For 
RX J0852.0-4622 Enomoto et al. (2006) and Aharonian et al. 
(2007) did not rule out either mechanism. 

The only well understood non-thermal emission mecha- 
nism is currently Synchrotron radiation of accelerated elec- 
trons. Synchrotron radiation is primarily observed in radio 
frequencies and is observed in X-rays in a few known SNRs 
(see e.g. Bamba et al. 2005, and refferences within). Non- 
thermal radio and X-ray observations are crucial for studying 
the non-thermal electron population. In the few known exam- 
ples of SNRs emitting X-ray Synchrotron radiation, the X-ray 
flux (per logarithmic frequency) is decreasing, indicating that 
the flux peaks at lower, unresolved photon energies. The radio 
and X-ray luminosities, and the implied position of the cutoff 
in the X-ray spectrum, may be used (and have been used in 
the models discussed above) for constraining the accelerated 
electron distribution, and thus for constraining the expected 
IC emission. 

In this paper we derive simple analytic relations between 
the dominant radio, X-ray and 7-ray emission mechanisms 
that can be used for distinguishing between IC and PP origins 
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of the 7-rays in SNRs were 7-rays were observed, and for 
predicting the 7-ray and neutrino flux values for SNRs where 
only radio and/or X-ray emissions where detected. When- 
ever possible, the simple analytic approximations we find are 
compared to, and shown to agree with, previous detailed cal- 
culations, with the advantage of being easier to follow and 
maintaining the explicit dependence on the unknown parame- 
ters. 

First we compare in section § 2 the expected IC, PP, Syn- 
chrotron and Thermal Bremsstrahlung (TB) fluxes in a simple 
one zone model of shocked ISM plasma. The plasma is as- 
sumed to consist of thermal and accelerated electron and pro- 
ton components with the later consisting of relativistic parti- 
cles having a power law distribution in energy. In this section, 
the high energy cutoffs of the accelerated electron and pro- 
ton distributions are ignored. Next, we discuss in section § 3 
the maximal energy attainable by electrons and protons in a 
SNR due to cooling and limited SNR age, and the implica- 
tions for the non-thermal emitted spectra, assuming diffusive 
shock acceleration (DSA) as the acceleration mechanism. We 
then find in section § 4 an upper limit to the value of K ep , the 
ratio of the number of accelerated electrons and the number of 
accelerated protons at a given energy, K ep > 10~ 3 by studying 
the radio observations of SNRs in M33. This parameter enters 
into the ratios of IC and Synchrotron emission to PP emission. 
This lower limit is used to rule out previously suggested SNR 
broadband emission models that used considerably lower val- 
ues. In section § 5 we apply the results of earlier sections to 
show that the broad-band spectrum of the SNRs RX J 17 13.7- 
3946 and RX J0852.0-4622 is inconsistent with a PP origin 
and is consistent with an IC origin of the > 1 TeV emission. 
We discuss previous claims that this emission cannot be due 
to IC and argue against them. The results are summarized and 
discussed in § 6. 

2. PP EMISSION VS. THERMAL AND NON THERMAL 
ELECTRONIC EMISSION 

In this section we compare thermal and non thermal con- 
tinuum emission mechanisms in the shocked plasma behind 
SNR blastwaves. The non thermal emission is assumed to 
be emitted by relativistic, accelerated electrons and protons 
with power law distributions in energy. In this section we 
ignore the energy cutoffs of the accelerated particle distribu- 
tions. This issue is discussed in § 3. We focus on ratios of 
the expected fluxes which are weakly dependent on unknown 
parameters such as distance to the remnant and total energy. 

First we write down in § 2.1 simple expressions for the lu- 
minosities due to the different processes in simple forms that 
allow easy comparison with each other (a derivation of these 
equations is given in § A). Next, we compare in § 2.2 the two 
~ 1 TeV 7-ray emission mechanisms, IC and PP. We then de- 
rive in § 2.3 constraints on 7-ray PP and IC emission by com- 
paring them to thermal X-ray Bremsstrahlung and to radio 
synchrotron emission. Finally, the results of § 2.1-§ 2.3 are 
compared in § 2.4 to earlier studies of the SNRs RX J 17 13.7- 
3946 and RX J0852.0-4622 (the only shell-type SNRs that are 
known to emit > 1 TeV 7-rays). 

We note that most of the results presented in this section are 
not restricted to SNRs and are applicable to any system that 
efficiently accelerates protons and/or electrons to relativistic 
energies with power-law energy distributions. 

2.1. Emission mechanisms 



Consider the shocked plasma in the downstream of the 
blastwave of a SNR. Here we consider radiation emitted by 
four distinct particle populations: 

1. Thermal electron and proton components with simi- 
lar number densities, n e <~ n p = n, which we assume 
consist of most of the particles. For simplicity we as- 
sume that the electron and proton energy distributions 
are given by Maxwelians with temperatures T e and T p 
respectively with T e = CeT p . The total number of pro- 
tons or electrons is N and the total thermal energy is 
E th « (3/2)NT p . 

2. Power law distributions of relativistic accelerated elec- 
trons and protons with an electron:proton ratio K ep , 
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The distribution of the protons is described by (1) for 
proton energies e p . m i n <e p < s p . max with e p „ in — m p c 2 . 
The value of e max depends on the SNR parameters and 
acceleration mechanism. Estimates of e max assuming 
DSA will be derived in § 3. 

We study the following radiation emission mechanisms: 

1. 7-rays and neutrinos emitted as a result of proton- 
proton collisions (PP) between the relativistic protons 
and the thermal protons. The PP gamma-ray luminos- 
ity per logarithmic photon energy is given by [cf. Eq. 
(A8)]: 

dN„ 



vh v p P = C w {p)2e p (v)'-^-a 1 p fnchv 
ds p 



(3) 



where e p dN p /de p is to be evaluated at e Pi pp(u) = lOhis 
(photon energies are referred to through the photon fre- 
quency throughout the paper), the typical proton en- 
ergy for which photons with energy hv are emitted. For 
p = 2,2.2 we have C PP (2) w 0.85,C PP (2.2) w 0.66. The 
neutrino luminosity is similar to the 7-ray luminosity at 
equal photon and neutrino energies. 

2. 7-rays emitted by Inverse Compton (IC) resulting from 
the interaction of the relativistic electrons with CMB 
photons. The IC gamma-ray luminosity per logarithmic 
photon energy is given by [cf. Eq. (A 12)]: 

1 dN e 4 2 

vL v ic = Cic(p)-e e — -<7t%(v)Ucmbc, (4) 

2 de e 3 

where 7cmb > Uqmb = fl 7c MB are the temperature and en- 
ergy density of the CMB photons. e e dN e /de e is to be 
evaluated at e e jc(v) = 7 e (^)m ( ,c 2 = m e c 2 (hv /3Tcmb) 1 ' 2 , 
the typical electron energy for which electrons up scat- 
ter CMB photons to energy hv. The correction factor, 
Cic(p), is approximately C\ C {p) w 0.8 (to within 5%) 
for 2 < p < 2.2. It is useful to note that j 2 (v)Ucmb = 
[^cmb/(37cmb)]^ ~ 0.9n C MB^ where h C mb is the 
number density of CMB photons. 
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3. Radio and X-ray synchrotron (Syn) emission of the 
relativistic electrons in an assumed magnetic field B. 
The synchrotron luminosity per logarithmic frequency 
is given by [cf. Eq. (A18)]: 



and 



vl}j* = C^{p)^e e {u)-^-a T l 2 e {v)U B c 



(5) 



where Ub = B 2 /(8ir). e e dN e /de e is to be evaluated at 
£e( v ) = "f e (v)m e c 2 = (2i>/i>B) l ' 2 m e c 2 , the typical energy 
of electrons emitting photons with frequency u, where 
vb = qB/(2Trm e c). The correction factor, Cs yn (/?), is ap- 
proximately Cs yn (p) s=s 0.8 (to within 5%) for 2 < p < 
2.2. 

4. Thermal-Bremsstrahlung (TB) emission of the thermal 
electrons interacting with the thermal protons. The 
maximal TB luminosity per logarithmic frequency is 
emitted at the photon energy hv = T e and is given by 
[cf. Eq. (A3)]: 

vL l B hv=T e = \J~^e~ l a e g ff N(T T nc^ m e c 2 T e (6) 

where e is the natural logarithm, a e w 1/137 is the fine 
structure constant and gjj is the thermal Gaunt factor. 
For 100 eV < T e < 10 keV, the value off// (for hv = T e ) 
is in the range, 0.8 < g/f < 1.2 (e.g. Karzas & Latter 
1961). 

We note that the amount of secondary electrons and 
positrons resulting from PP interactions is most likely neg- 
ligible compared to the primary population of accelerated 
electrons. The energy output in electrons and positrons 
per logarithmic particle energies is roughly equal to the 
7-ray emission given by Eq. (3). The ratio of sec- 
ondary electrons+positrons to protons for an SNR of age 
t = lOOOf^yj-yr evolving into a medium with proton den- 
sity n = iiq cm -3 is thus roughly given by (ignoring cool- 
ing, which affects both primary and secondary populations 
in the same way) e 2 dN e+e -/de ~ 0.2e 2 dN p /de p a p n p l nct <~ 
10~ 6 e 2 dN p / ds p t ^yj-rio. As long as K ep » 10~ 6 ? ky r «o> the 
contribution of the secondary electrons to the broad band 
emission is negligible. Henceforth we ignore this contribu- 
tion. 

2.2. IC to PP emission ratio 

Here we directly compare the two competing TeV 7-ray 
emission mechanisms. Ignoring the possible cutoffs of the 
spectrum of both species, the ratio of expected IC to PP 
gamma ray luminosities per photon frequency can be approx- 
imated by [compare Eq. (3) and Eq. (4)]: 
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where K ep =10 2 K ep -2 and n = l«o cm 3 . It is useful to note 
that 
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Comparison of Eq. (7) with the results of previous studies is 
presented in § 2.4. 

An electron to proton ratio of order K ep <~ 10" 2 is commonly 
assumed based on the measured electron:proton ratio in the 
cosmic rays (see e.g. Longair 1994) under the assumption 
that SNRs are the main source of proton and electron cosmic 
rays. In section § 4 we find a lower limit of K ep > 10~ 3 based 
on radio observations of SNRs in M33. 

Using Eq. (7), we see that as long as electron cooling does 
not suppress the IC flux, IC dominates PP emission as long as 



n< 10^_ 2 4y )/2 cm- 3 . 
The effect of electron cooling is addressed in § 3. 
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2.3. PP and IC to TB X-rays and non thermal radio Syn 

Here we compare the expected PP and IC 7-ray emission 
to thermal and synchrotron emission. This is useful for con- 
straining the expected gamma-ray and neutrino fluxes based 
on observed radio and X-ray fluxes. 

By comparing equation (3) to (6) we see that the ratio of PP 
7-ray luminosity at photon energies hv 1 to the X-ray TB lumi- 
nosity at photon energies hv x = T e (the photon energy of max- 
imal emission per logarithmic photon energy) can be written 
as: 

v~,L, 
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where w (3/ 2)NT p is the total thermal energy. Using equa- 
tion (1) this can be written as: 
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where e p = 0.U Pr i = E p /E th , A p = l0A Ptl , T p = T pM y keV, 



hv 1 = i>iev TeV and we substituted e„ 



Assuming 



p > 2, e max > 10 TeV and hv^ > 10 GeV the factor in the 
second line of Eq. (12) is smaller than 1.1. 

Temperatures T p > keV are characteristic of the shocked 
plasma in young SNRs with blastwaves propagating at veloc- 
ities v s > 1000 km s -1 . In fact, the proton temperature behind 
a strong shock propagating with velocity v s is given by: 



T 
1 p 



■■ 2vi keV, 



(13) 



where v s = 1000v 8 km s" 1 and an adiabatic index equal to 7 = 
5/3 was assumed. A lower limit to the shock velocity and 
thus to T p for SNRs where non-thermal X-rays are observed 
is discussed in § 3. 

The ratio ( e of electron to proton temperatures depends 
on the amount of collisionless heating in the shock and the 
following heating through Coulomb scattering. The amount 
of collissionless heating for high Mach shocks is not really 
known (for a recent review see Rakowski 2005). A lower 
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limit to ( e can be derived by assuming that there is no colli- 
sionless heating. After a time t = fk yr kyr the ratio would be 
[cf. Eq. (A5)]: 



> 0.6 (A e/ ,.i. 5 «ofkyr) 2,/5 T^feV' 



(14) 



where X ep = 30\ epA 5 is the Coulomb logarithm. Eq. (14) was 
derived assuming m e /m p <( e <l and is valid as long as the 
resultant value is in this range. Substituting Eq. (14) in Eq. 
(12) we find: 
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By comparing equations (3) and (5) we see that the ex- 
pected ratio of 7-ray PP luminosity at photon energies hv^ 
to the radio synchrotron luminosity at frequency u R can be 
approximated by: 
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where 7 e ,s yn (^K) = e e ,s y n(^)/ OT e c = (A-Km e cv R / qBy' 2 is the 
typical gamma factor of electrons emitting radiation with 
vr = VQYh. GHz frequency and B = IOB-5 fiG. Assuming p > 
2, the factor in the second line of Eq. (16) is smaller or equal 
to 1 . It is useful to note that, 



= 10 ^ h ^ 2 „ 2 X lOVf^VeV (17) 
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Comparison of Eq. (16) with the results of previous studies is 
presented in § 2.4. 

Ignoring the possible cutoff of the IC spectrum, the ex- 
pected ratio of TeV IC emission to GHz synchrotron emis- 
sion is approximately given by [compare Eq. (4) and Eq. (5)]: 
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The luminosity ratios of the different radio, X-ray and 7- 
ray emission mechanisms are given by equations (7),(12), 
(16) and (18). Flux normalization is obtained by noting that 
the expected radio flux per logarithmic frequency for a SNR 
with total energy E = 10 51 £5i erg and a fraction rj p = 0.1ry p _i 
of the total energy carried by accelerated protons E p = rj p E 
(f] p <~ e p /2) located at a distance d = d\^ c kpc, is approxi- 
mately given by [using Eq. (5)], 



^/iASyn|iy=lGHz ' 
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Comparison with previous studies 

Next we compare the results presented in this section to 
previous studies of the broad-band emission of SNRs. We 
focus on studies that were published following the discovery 



of the > 1 TeV 7-rays from the shell-type SNRs RX J1713.7- 
3946 and RX J0852.0-4622 (the broad-band emission of these 
SNRs is discussed in § 5). Recent broad-band studies of RX 
J1713.7-3946 where done in (Aharonian et al. 2006; Berezhko 
& Volk 2006; Porter et al. 2006; Moraitis & Mastichiadis 
2006) while studies of RX J0852.0-4622 include (Enomoto 
et al. 2006; Aharonian et al. 2007). These studies differ in 
the way the particle distributions are obtained, in the assump- 
tions regarding the magnetic field value and in the assump- 
tions regarding the ambient IR radiation field. Berezhko & 
Volk (2006) numerically solved time-dependent CR trans- 
port equations, coupled nonlinearly with the hydrodynamic 
equations for the thermal component. Aharonian et al. (2006), 
Porter et al. (2006) and Aharonian et al. (2007) assumed a 
constant injection of particles with a power-law spectrum that 
is cutoff exponentially over a fixed period of time. They cal- 
culated numerically the effects of cooling on the particle spec- 
trum (Aharonian et al. 2007, take in addition particle escape 
into consideration). Moraitis & Mastichiadis (2006) found an 
analytic solution to "two-zone" (acceleration zone and escape 
zone) spatially averaged kinetic equations that include cool- 
ing. Enomoto et al. (2006) assume a power-law spectrum that 
is cutoff exponentially. Berezhko & Volk (2006) estimated 
the value of the magnetic field based on observation of thin 
X-ray filaments (see discussion in § 5.4. 1) while the other au- 
thors allowed for different magnetic field values. Porter et 
al. (2006) included a detailed model of the galactic radiation 
field with IR and CMB dominating in different places, while 
the other authors assumed 'standard' averaged values. 

All the above studies focused on the non-thermal emission 
mechanisms only. Comparisons of Eqs. (7) and (16) with the 
results of these studies are presented in tables 1 and 2 respec- 
tively. As can be seen, there is good agreement (up to a 



TABLE 1 

Ratio of IC to PP emission for negligible suppression of IC 

DUE TO COOLING 



Ref. 
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n[ cm 3 ] 


L v ic/K PP ( TeV) 


Eq. (7) 
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The value of L„ \c/L v pp( TeV) that results from Eq. (7) for the values of 
K ep and n[ crrT 3 ] used in each reference is shown in the last column next to 
that obtained in each reference (fourth column). References: (1) Moraitis & 
Mastichiadis (2006); (2) Aharonian et al. (2007), figure 17a; (3) Aharonian 
et al. (2007), figure 17b; (3) Enomoto et al. (2006); 

The line distinguishes between PP and IC dominated models of SNRs RX 
J1713.7-3946 and RX J0852.0-4622. 

a The value of K ep was calculated by K ep S3 2oPq 2 'fio'Po*' 2 assumm g 
the escape zone dominates and that p, po <C p ma x, using the authors notations. 

factor <~ 2) between our analytic expressions, Eqs. (7) and 
(16), and the results of earlier detailed numerical calculations 
of the remnants RX J1713.7-3946 and RX J0852.0-4622. 

3. ENERGY CUTOFFS 

In section § 2 we discussed the radiation emitted by power- 
law distributed electrons and protons. In reality, the particle 
distribution functions can be approximated by a power law 
function only over a limited range of particle energies. The 
maximal energies and corresponding cutoff frequencies in the 
emitted spectrum were extensively studied before (see e.g. 



TABLE 2 

Ratio of PP to synchrotron emission 
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The value of vL v pp( TeV)/i/L„ s yn ( GHz) that results from Eq. (16) for the 
values of K ep , n[ cm -3 ], B[ /iG] and p used in each reference is shown in 
the last column next to that obtained in each reference (sixth column). Ref- 
erences: (1) Berezhko & Volk (2006); (2) Moraitis & Mastichiadis (2006); 
(3a) Aharonian et al. (2007), figure 18a ;(3b) Aharonian et al. (2007), figure 
18b; (3c) Aharonian et al. (2007), figure 17a; (3d) Aharonian et al. (2007), 
figure 17b; (4) Enomoto et al. (2006) 

The line distinguishes between PP and IC dominated models of SNRs RX 
J1713.7-3946 and RX J0852.0-4622. 
VL^pptTeV)/^ syn(GHz) 

b In this reference the ambient density is nonuniform, n is taken as the 
value of the ambient number density currently encountered by the shock 
c See footnote a in table 1 

Druryetal. 1994; Reynolds 1998). For completeness we write 
down in this section the expressions for the maximal particle 
energies attainable by DSA in a simple SNR model and the 
possible spectral cooling break in the electron spectrum and 
discuss the implications for the spectrum of the emitted radi- 
ation. 

First we write down in § 3. 1 the maximal energies attainable 
by DSA as a function of the SNR radius, age and energy, ig- 
noring the dynamical relation between these quantities. Next, 
we focus in § 3.2 on the Sedov-Taylor (ST) SNR evolution 
phase. We then discuss in § 3.3 SNRs in which non-thermal 
X-rays are observed. We find a lower limit to the shock ve- 
locity and post shock temperature in such SNRs. In addition 
we derive constraints on the 7-ray and X-ray spectral cutoffs 
that must be satisfied by an IC model for the 7-rays emit- 
ted by such SNRs. Finally, we find in § 3.4 an upper limit 
to the PP emission for SNRs in which the IC TeV emission 
is suppressed by synchrotron cooling of the energetic elec- 
trons. This is done by comparing the PP emission to the radio 
synchrotron emission with the implied minimal value of the 
magnetic field that is required to cool the electrons in times 
shorter than the SNR age t. 

3.1. Energy cutoffs 

Consider a SNR with the following parameters: Energy in 
shocked matter E = 1 5 1 £5 1 erg, radius R = 10R 1 pc, age t = 
fky r kyr shock velocity v, = lOOOvg km s~\ ambient medium 
density of n = «o cm" 3 at a distance of d = d^ pc kpc. The dis- 
tance to the SNR is related to the radius by d^ pc « R\/0° 
where 9° is the angular diameter of the SNR on the sky in 
degrees. The shock velocity, age and radius are related by 
\ s = aR/t w \(faRit^ r cm s" 1 with 0.4 < a < 1, the lower 
limit obtained for Sedov-Taylor (ST) expansion and the upper 
limit for free expansion (FE). 

We assume that electrons and protons are accelerated to 
power law spectra [cf. Eq. (1)] with p w 2 up to cutoff en- 
ergies, £ e . max and e^. max respectively, with a possible cooling 
break in the electron spectrum at m p c 2 < £ e ,break < £«,max be- 
yond which the power law index is p + 1 . 



Both electron and proton energies are limited by the finite 
available acceleration time due to the finite SNR age. The 
maximal proton or electron energy due to the finite time sat- 
isfies: 

f a cc(£max) = t, (20) 

where t acc (e) is the time it takes electrons or protons to reach 
energy e (assumed to be equal for protons and electrons). The 
maximal energy of accelerated electrons can also be limited 
by cooling, in which case we have 



^acc(^^.max) — ^cool(^e.max); 



(21) 



where f C ooi(£e) is the cooling time of electrons with energy 
e e . In the later case, a cooling break is expected at an energy 
£ e ,break satisfying 

fcoolfe.break) = t . (22) 



The acceleration time, t- icc (e), can be approximated by, 

^ ^cycle 



5 A/3' 



(23) 



where A/3 = (v s ~ u d)/ c < u d ~ v s/4 is the downstream veloc- 
ity and f cyc ie is the shock crossing cycle time. Assuming the 
downstream residence time dominates the cycle time, we have 



^cycle — 



4Dd 

UdC ' 



(24) 



where Dd is the downstream diffusion coefficient. The diffu- 
sion coefficient can be expressed as 



3qB 



(25) 



where £ is a dimensionless coefficient that satisfies £ > 1 with 
£, = 1 obtained in the Bhom diffusion limit. Using equations 
(23)-(25), we can write the acceleration time as 



16 e 
T qBf3 s w s ' 



(26) 



Using Eqs. (20) and (26), the maximal energy due do to the 
limited SNR age can be expressed as 

e^fecc =t) w ~ 60 r'aft-svg*! TeV. (27) 

The electron cooling time due to synchrotron emission is 
given by 



fcool — 



(4/3)^a T cU B 



HwkvY 1 ** 1 * (28) 



[we neglect the effect of IC which corresponds to a cooling 
time of ss 100(e«/10 TeV)" 1 kyr which we assume is much 
larger than the SNR age]. We thus expect a cutoff to the syn- 
chrotron and IC spectra at photon energies given by 

^Synfecc = ~ lB^V^CV keV, 

/^icfecc = t) ~ lOfiVs^CV TeV, (29) 
and a cutoff to the PP spectrum at photon energies given by 

hvMtzcc = t) ~ lOfi-sVsfliC 1 " TeV. (30) 

Electron cooling will be relevant at the energy given by Eq. 
(27) for strong enough magnetic fields, 



B> 10 



a 2 R\ 



1/3 



/xG. 



(31) 



If Eq. (31) is satisfied, there would be a spectral break at syn- 
chrotron and IC photon energies given by: 

hvs yn {t mo \ =t)fts 3Bl5% y 2 r keV, 

^cfcooi = t) « 40BI^ k - y 2 r TeV, (32) 

and a cutoff at photon energies of 

3 3 

^Synfcool = face) ~ l Ot e l m e V* ~ £ l 0.15v£ keV, 

hi^citcooi = face) ~ C^vfflj TeV. (33) 

The synchrotron and IC flux per logarithmic photon energy 
at photon energies above the break and below the cutoff would 
be suppressed by a factor of 



vL v (with cooling break) 
vL v (no break) 



hv 



hv{t coo \=t) 



-1/2 



(34) 



compared to the flux that would be emitted by a power-law 
without a break. 

3.2. SNR dynamics 

We now focus on the Sedov-Taylor (ST) phase. The SNR 
enters the ST phase when the mass of the swept up ambient 
medium, 



M swep t~ lOO^n o M , 
is larger than the ejecta's mass. In this case we have: 
a = 0.4, 



(35) 



0.5pR 5 



-Sflfnf^kvr, 



lO^VV^ms" 1 . 



E 

v s = OAR/t < 
Substituting Eq. (36) in Eqs. (29)-(33) we obtain: 

^Synfcool = OsT ~ 0. 15E 5 iBZsR{ 5 Hq 1 k eV, 
^Syn(f C ool = facc)sT ~ O.lf^si^/lQ 1 keV, 
/^Synfecc = OsT ~ 0. 1^-5^5 1*T V keV, 

fc^cfcooi = Ost ~ 2B:^5i^T 5 «o 1 TeV, 

hlAdtcool = facc) ST - 1.5^^51*1 V TeV, 



(36) 



hv\cih< 



:f)sT-ir 2 B-5£5i«T 1 no 1 TeV, 



and 



/^Ppfecc = OST ~ 3^-3^51 



1 /2 J? -l/2„-l/2 



TeV. 



(37) 

(38) 
(39) 



3.3. SNRs with observable non-thermal X-rays 

For SNRs with observable non-thermal synchrotron X-rays, 
we can find a lower limit to the shock velocity by demanding 
that there will be no cooling cutoff for photons with energies 
smaller than hv x = iA t w keV, i.e. hvs yn (t coo \ = f acc ) > hv x . Us- 
ing Eq.(33), this can be written as 



v s > 3 x 10 8 £ 



1/2 1/2 



cm s 



(40) 



The minimal velocity constraint has several implications. 
First, this can be used to obtain a minimal value for the pro- 
ton temperature in the downstream. Comparing Eqs. (33) and 
(13), we find: 



t p > „£ a e — hi>x ~ I0£,hv x . 
9 m e 



(41) 



Second, assuming that the shock velocity is not much larger 
than 3000V8.5 km s" 1 , the diffusion coefficient cannot be 
much larger than the Bohm limit (£ = 1), 



(42) 



Under this assumption, the proton temperature is constrained 
by: 

lOZhvx <T P < 20v| 5 keV. (43) 

Third, using E > 3pv^R 3 (which is valid for both the ST 
and FE phases), we find a lower limit to the ambient medium 
density of 



" <0 - 1 ^^keV Cm " 3 - 



(44) 



Next we compare the cutoff in the IC emission to the cutoff 
in the synchrotron radiation. The energies of photons emit- 
ted by electrons through IC and Synchrotron are both pro- 
portional to the square of the Lorentz factor of the emitting 
electrons. The ratio of photon energies emitted through IC by 
electrons to the photon energies emitted through synchrotron 
by the same electrons is approximately given by: 

-T^~3T C mb ^-«10 10 flj. (45) 

The ratio of IC to synchrotron power, emitted by the same 
electrons, is approximately given by: 



0.1BZ 



(46) 



In particular, the photon energies where the IC and the syn- 
chrotron luminosities are cutoff should satisfy Eq. (45), and 
the luminosity values at these photon energies should satisfy 
Eq. (46) (this is true in principle for any feature in the spec- 
trum). We can use both equations to write a constraint that 
does not depend on the value of the magnetic field (or the ac- 
celeration mechanism): 



^/C J cutof L ^3 xl()1 o 



'vie I ^/c.cuioff 



hv 



'Syn, cutoff 



(47) 



A note of caution is in order regarding equation (47). A 
'cutoff frequency is not a well defined quantity in general. 
For known functional forms, prescriptions for defining a spe- 
cific frequency can be given. The precise value of the numeri- 
cal coefficient in (47) may be somewhat different for different 
prescriptions. In addition it should be noted that while the 
IC spectrum of a single electron has a sharp cutoff (photons 
with energies larger than the initial electron energy cannot be 
generated), the synchrotron spectrum cuts-off exponentially, 
resulting in different photon spectra for given cutoff forms. 
Taking this into consideration and since the precise electron 
spectrum is not known, the cutoff frequencies are defined only 
to within an order of magnitude. 

3.4. Suppression ofIC due to radiative cooling 

The calculated expected ratio given by Eq. (7) is valid 
as long as there is no significant suppression of the elec- 
tron population due to cooling. Electrons responsible for 
TeV IC emission have Lorentz factors of approximately 

7 ~ ( TeV^ e v/3rcMB)^ 2 ~4x lO 7 ^]/ 2 , an d a corresponding 
cooling time of [cf. Eq.(28)]: 



fcool 



= 6v T % 2 B_ 2 5 kyr. 



(48) 
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Cooling will affect these electrons only if the cooling time 
is shorter than the lifetime t of the SNR which would be true 
only if the typical magnetic field is large enough: 



1/2 r 
TeV 'kyr 



(49) 



A larger magnetic field would imply stronger synchrotron 
emission. We can use this to write a constraint on the TeV 
PP emission in case the IC emission is suppressed. Assuming 
that the electrons responsible for the IC TeV emission were 
suppressed by cooling, we can use Eqs. (49) and (16) to ob- 
tain: 



VfL v ^ pp 

VrL ur Syn 



(50) 



Another constraint can be derived by comparing the PP 
emission to the Synchrotron radiation at X-ray frequencies 
assuming that the electrons emitting the X-rays are also af- 
fected by cooling. This assumption is reasonable since we 
assume that electrons responsible for TeV IC emission are af- 
fected by cooling, and these electrons are responsible for syn- 
chrotron radiation of photons with energies hv > 100fi_5 eV 
[cf. Eq. (45)]. Using Eqs. (16) and (34), the ratio of the 
PP flux to the X-ray flux in the frequency range between the 
cooling break given by Eq. (32) and the cooling cutoff given 
by Eq. (33) (for a spectrum with p > 2 this is the maximum 
value of vLvsy n ) is given by 



v^L Vi pp 

v xL Vx Syn, max 



■ 0.4« a la fc tK~ l ~ 1.5 x KT^^kyrno, 



(51) 

where we assumed p = 2. This equation has a weak depen- 
dence on p since the X-ray emitting electrons have energies 
that are similar to the TeV 7-ray emitting protons. Equation 
(51) has the following simple interpretation. Suppose that the 
amount of protons per unit energy and unit time that are be- 
ing accelerated by the shock is given by Q(e). The amount of 
protons per unit energy at an age t is roughly dN/de <~ Q(e)t 
and so the PP luminosity per logarithmic frequency is roughly 
given by [cf. Eq. (3)] vLv ?? ~ 0.2Q(e)s 2 tnca™\ The 
electron injection rate at electron energies of e is K ep Q(e). 
As the electrons are constantly being cooled, the energy in- 
put in accelerated electrons is equal to the energy emitted 
in synchrotron radiation. The X-ray synchrotron luminos- 
ity per logarithmic frequency is thus roughly vLi/syn coo i e d <~ 
0.5K ep Q(e)e 2 (the factor of 0.5 comes from the fact that the 
logarithmic interval in photon energies is twice that of the 
emitting electrons due to the v oc j 2 dependence). The ratio 
of these expressions is equal to the result in Eq. (51). 

4. LOWER LIMIT ON K EP FROM EXTRAGALACTIC 
SNRS 

In this section we find a lower limit for K ep using the ob- 
served radio fluxes from large SNRs in M33 assuming that 
the value of K ep does not vary significantly between SNRs. 

One way to to estimate the amount of accelerated electrons 
is through the radio synchrotron emission. The radio lumi- 
nosity is determined by the energy in accelerated electrons 
and by the magnetic field value. The amount of energy in 
accelerated electrons cannot be deduced if the value of the 
magnetic field is not known. An upper limit to the magnetic 
field is given by the requirement that the magnetic field does 



not exceed equipartition. Here we assume that the fraction eg 
of the thermal energy carried by the magnetic field behind the 
shock does not significantly exceed e# <~ 0.1. In addition we 
assume that the fraction r} p ~ e p /2 of the total energy carried 
by relativistic protons does not significantly exceed r\ p <~ 0.1. 
Using equations (1) and (5) we can approximate the expected 
luminosity at 1 GHz by: 



,3/2 



L,, S yn( GHZ) « 4 X 10 Z %,,_2%-l£ 5 lB^ 

x A-\(5B~_l /2 r (p - 2) erg cm" 2 s" 1 Hz" 1 (52) 

where i] p = 0.1 r\ p .-\- For an assumed maximal Lorentz factor 
7/j,max ~ 10 5 , the factor in the second line of Eq. (52) equals 
s=s 0.9 for p = 2 and w 1.6 for p = 2.2, and will be ignored 
henceforth. 

The ratio of the magnetic to thermal energies behind the 
shock in the Sedov-Taylor phase can be approximated by: 



B 2 /U 



B 



, ~ 3—R 1 E~ l «4x 10- 4 B 2 _ 5 R]E^ . (53) 

o7T 



Extracting the magnetic field from Eq. (53) and substituting 
it in Eq. (52) we have 



erg s 1 Hz 1 . 
(54) 



L^ Syn ( GHz) « 3 x 10 2 X,-2^,-i£5i /4 4^i^ 9/4 

Radio luminosities of SNRs with known distances in nearby 
galaxies (including the milky way) are summarized by Ar- 
butina et al. (2004), and virtually all have luminosities greatly 

exceeding 3 x \Q 22 R^ A erg s" 1 Hz -1 , the typical value ex- 
pected from Eq. (54) for K ep = 10~ 4 . However, we should 
stress that it is dangerous to reach conclusions based on such 
comparisons, since the observed luminosities are limited from 
below by the detectors' sensativities. Here we focus on a sam- 
ple of SNRs in M33 which is perhaps the most complete sam- 
ple of radio SNRs with known distances in a single galaxy 
(Gordon etal. 1999). 
Using Eq. (35) we see that SNRs with radii larger than 



R>2 



( ggj 

ViOM 



1/3 



-1/3 

P c 



(55) 



are in the ST phase. The smallest SNR in the sample has a 
radius of R ~ 5 pc and most SNRs in the sample have radii 
R > 10 pc. It is thus reasonable to assume that the SNRs in 
the sample are in the ST expansion phase. In fact, Gordon et 
al. (1998) have shown that the radii distribution function of a 
larger optical SNR sample that includes the radio SNR sample 
is consistent with ST expansion (and is inconsistent with free 
expansion). 

The luminosities of the observed SNRs in M33 are shown in 
figure 1 along with the observational threshold (dashed line) 
and the expected limits according to Eq. (54) correspond- 
ing to K ep = 10" 4 (lower, green) and K ep = 10" 3 (higher, red), 
adopting a distance of d = 840 kpc to M33. 

As stressed by Gordon et al. (1999), there are probably un- 
observed SNRs with luminosities that fall beneath the obser- 
vational threshold. In fact, there are about twice as many 
SNRs seen in optical wavelengths (Gordon et al. 1998, the 
factor being roughly radius independent, e.g. 8, 23 and 34 
SNRs in the radio sample with radii R < 10, 15, 20 pc respec- 
tively compared to 15,42 and 67 SNRs respectively in the op- 
tical sample). Assuming the optical sample is not far from 
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FIG. 1.— Radio 20cm luminosity of SNRs in M33 (Gordon et al. 1999). 
The dashed line is the observation flux (total and density) limit. The full lines 



are given by Eq. (54) with K ep - 
red). 



10 (lower, green) and K el 



■■ 10" 3 (higher, 



completeness, it is reasonable that roughly half of the SNRs 
are missed in the radio sample (this is true for R < 20 pc, 
while for R > 20 pc the optical sample is probably incom- 
plete, Gordon et al. 1998). Still, it is quite clear from figure 
1 that the luminosity implied by a value K ep = 10~ 4 is lower 
than the typical luminosity of large remnants by at least an 
order of magnitude. As an illustration, in order to reconcile a 
10~ 4 with the four SNRs observed with radius 

-l 



value of K, 



ep 



R ~ 50 pc and luminosity L u ( GHz) ~ 10 erg s Hz , their 
energies would have to be unreasonably high: 



E ~ 5 x 10 3 



R 



50 pc 



9/7 



K 



-4/7 ~V7 -3/7 
-4'lp-l e B,-l 



ep 



L v ( GHz) 



10 24 erg s" 1 Hz" 



4/7 



erg. 



(56) 



Note that SNe with energies that are larger than 10 52 erg 
(termed Hypernovae) have been detected (see e.g. Nomoto et 
al. 2006, and references therin). However, the estimated frac- 
tion of core-collapse SNe that belong to this group is of order 
10" 3 (Podsiadlowski et al. 2004) and thus having 4 SNRs with 
energies exceeding 10 52 erg among the ~ 100 SNRs in M33 
is unlikely. 

Using Eq. (54), we conclude that K ep > 10~ 3 is a reasonable 
lower limit and that K ep ~ 10~ 4 can be conservatively ruled 
out. 

A possible caveat in the arguments in this section comes 
from the fact that it is possible that the ambient CR electrons 
that have been swept up by the shock have a considerable con- 
tribution to the Synchrotron emission (Anderson & Rudnick 
1993). The arguments in this section will nevertheless remain 
valid in this case too, provided the ratio of the accelerated 
electrons and protons populations, including the CR contribu- 
tions, are similar for different SNRs. We note that assuming 
that the cosmic rays that enter the shock are reaccelerated by 
DSA, the shape of the spectrum of the population of relativis- 
tic particles will approach a power law and will not be affected 
by the distribution of the CRs in the ISM (Drury et al. 1994). 

5. APPLICATION TO RX J 17 13. 7-3946 AND RX 
J0852.0-4622 

Here we apply the results of sections § 2-§ 4 to show that the 
broad-band spectrum of the SNRs, RX J1713.7-3946 and RX 
J0852. 0-4622 is inconsistent with a PP origin and is consistent 



with an IC origin of the > TeV emission. First we summarize 
in § 5.1 the broadband observations of these SNRs. Next, we 
show in § 5.2 that a PP source for the observed 7-ray flux is 
inconsistent with the broad-band emission in these SNRs. We 
then show in § 5.3 that an IC source for the observed 7-ray 
flux is consistent with all observations. We show that the con- 
tribution of the PP 7-ray emission is negligible and argue that 
the neutrino emission from these SNRs is probably too low to 
be detected by current and planned neutrino telescopes. Fi- 
nally, we compare in § 5.4 the results presented here to previ- 
ous studies. In particular we discuss previous claims against 
an IC source of the 7-ray s. 

5.1. Characteristics ofRX J 17 13. 7 -3946 and RX 
J0852. 0-4622 

The observations of these SNRs are described by Aharonian 
et al. (2006, 2007) and references therein. Some of the main 
features are summarized below. In many ways these two shell 
type SNRs are similar . Both have comparable radio and TeV 
fluxes, 

vf v |ghz ~ few 10" 13 erg cm" 2 s" 1 , 



vfv\TeM ~ few 10 11 erg cm 2 s 



(57) 



span similar angles on the sky (6 1,2° respectively), have 
non thermal X-ray emission, which is consistent with a cut- 
off frequency of order hv cuto ff < keV. The gamma ray energy 
flux is consistent in both SNRs with a power law vf v oc v° and 
an exponential cutoff at photon energies of ~ 10 TeV (for RX 
J0852. 0-4622 the detection of the cutoff is less clear, Aharo- 
nian et al. 2007). 

Perhaps the main difference is in the hv ~ 1 keV X-ray flux 
which is larger for RX J17 13.7-3946 by a factor of about 5, 
vf v \ ps few 10~ I0 ,few 10~ n erg cm" 2 s" 1 respectively. 

One of the main characteristics of these SNRs is a non 
thermal dominated X-ray emission. The lack of observable 
thermal radiation can be used to obtain an upper bound on 
the value of the ambient density. Number densities consider- 
ably smaller than 1 cm" 3 where obtained (Slane et al. 1999, 
2001; Pannuti et al. 2003; Cassam-Chena'i et al. 2004), which 
in turn constrain the amount of proton-proton collisions. For 
RX J1713. 7-3946, limits on n from lack of thermal radiation 



of n < 0.3 cm-Vkpc/6)" 1 / 2 , n « 0.05-0.07 cm"V kp c/6)' 



-1/2 

and n < 0.02d k ^ 2 cm" 3 were obtained by Slane et al. (1999); 
Pannuti et al. (2003) and Cassam-Chenai et al. (2004) re- 
spectively. For RX J0852. 0-4622 a limit on n from lack 
of thermal radiation (for temperature greater than 1 keV) of 

n < 0.034 p ' c /2 cm" 3 was obtained by Slane et al. (2001). 

There have been claims that RX J 17 13.7-3946 is interacting 
with molecular clouds (Slane et al. 1999; Fukui et al. 2003, 
at 6kpc and lkpc respectively). Interaction with molecular 
clouds of both SNRs is unlikely given the low densities im- 
plied from lack of thermal radiation and the observed roughly 
homogenous emission (Aharonian et al. 2006, 2007). The 
positive TeV to CO line emission correlation that was claimed 
for RX J17 13.7-3946 is not convincing since the CO inten- 
sity changes by some two orders of magnitude while the TeV 
changes by a factor of 2 (average to peak, Aharonian et al. 
2006, 2007). In any case, interaction with molecular clouds 
cannot account for the entire emission and we will ignore this 
possibility henceforth. 

Distance and age estimates for these remnants are inconclu- 
sive (Aharonian et al. 2006, 2007, and references within). We 
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think that it is worth mentioning that claims that the distance 
to these SNRs is d < 1 kpc (e.g. Fukui et al. 2003; Aschen- 
bach et al. 1999) require some coincidence since the galactic 
latitude of both SNRs is b < 1°, (b = 0.5, 1.2 for RX J1713.7- 
3946 and RX J0852.0-4622 respectively) whereas the SNRs 
at this distance should be distributed in the range \b\ < 10° 
assuming SNRs are distributed homogenously throughout the 
galactic gaseous disk hight. For RX J 17 13.7-3946 the coin- 
cidence that is required is more extreme since this SNR lies 
in the direction of the galactic center, b = 0.5° , I = 347°, close 
to a 'hole' in the galactic CO line emission (Slane et al. 1999; 
Moriguchi et al. 2005). These positions on the sky may not be 
coincidental if these SNRs are farther away- a few kpcs from 
us (as most SNRs are). On the other hand, we note that such 
a coincidence is certainly possible and we do not assume in 
what follows that the distance to these remnants is larger than 
1 kpc. 

5.2. Upper bounds on PP emission 

We first consider the constraints on the PP emission result- 
ing from the comparison of the PP emission to the IC and 
Synchrotron non-thermal emission. By inserting K ep <~ 10~ 2 
and n < 0. 1 cm" 3 in Eq. (7) we see that unless the IC emission 
is suppressed by cooling, the > 1 Te V emission in these SNRs 
is completely dominated by IC. Irrespective of cooling, the ra- 
tio of PP > 1 TeV emission to the synchrotron radio emission 
is given by Eq. (16). A lower limit for the magnetic field is 
given by (49) for the case where IC emission is suppressed 
by cooling. Alternatively, a lower limit of B > 10 ^G can be 
derived by demanding that the IC emission generated by the 
electrons that emit the observed X-ray synchrotron emission 
does not exceed the observed gamma ray emission [using Eqs. 
(45) and (46), see e.g. Aharonian et al. (2006)]. By inserting 
n = 0.1n-i cm" 3 and B > 10B- 5 ^G in Eqs. (16) and (50) and 
assuming p>2,we find that 



V-yL^pp , 3/2 

^5K ep _ 2 B_ 5 n_i 



VGWzL VQYiL Syn 



(58) 



and 



V ^ PP < , K -l ,3/4 3/8 
fGHzi^ GHz Syn 



with the later equation applicable if the IC emission of 
photons with energy hv = TeV is suppressed by cool- 
ing. Comparing this to the observed ratio of fluxes 
per logarithmic frequency at hv = 1 TeV and v = 1 GHz, 
which for RX J1713.7-3946 and RX J0852.0-4622 is 
vf v ( TeV)/vf u ( GHz) ~ 100, we see that the contribution 
of the PP TeV emission is negligible compared to the total 
> 1 TeV emission. 

In case that the synchrotron X-rays are also affected by 
cooling (this is likely if the IC emission is suppressed by cool- 
ing), by using equation (51) we find that 



pp 



VxLi/ x Syn.max 



1.5 x 10 A K ep -2t^y r n-\ . 



(60) 



Comparing this to the observed ratio of fluxes per 
logarithmic frequency at hv = 1 TeV and hv <~ 1 keV, 
which for RX J17 13.7-3946 and RX J0852.0-4622 are 
vf v { TeV)/vf l/ ( keV) ~ 10 and 2 respectively, and assuming 
that the maximum Synchrotron luminosity cannot be much 



higher at lower frequencies, we see again that the contribu- 
tion of the PP TeV emission is negligible compared to the 
total > 1 TeV emission. 

We have assumed above thntK ep ~ 10" 2 , in accordance with 
the local ratio of CR electrons to protons and with section § 4. 
It should be emphasized that a direct estimation of the value 
of K ep using the arguments of section § 4 is not possible for 
these SNRs, since their radio luminosity is not known (due to 
the uncertain distances) and since the value of eg is not known 
for these remnants (the expected value of eg is discussed in 
§ 6). In § 5.3 it is shown using Eq. (19) that a magnetic field 
value of B ~ 10 /LtG, required in the IC scenario, and a value 
for K ep of K ep <~ 10~ 2 are consistent with the observed radio 
flux for a distance of <~ 1 kpc. 

We next consider the constraint on the PP emission result- 
ing from the comparison of the PP emission to the TB X-ray 
emission. The ratio of > 1 TeV PP Luminosity to TB Lumi- 
nosity is given by Eq. (12). 

Constraints on the shock velocity and more importantly the 
post-shock proton temperature for SNRs with observable non- 
thermal X-ray radiation are given by Eqs. (40) and (43) 



v s > 3 x 10 8 £ 1/2 ^ev cm s"\ 



and 



lO^ev < T p < 20vf. 5 keV 



(61) 
(62) 

respectively (the upper limit to T p results from the assump- 
tion v s < 3000v| 5 km s" 1 ), where we assumed that the cut- 
off in the X-ray spectrum is at hv cu i s = i\ e y keV. £ is the 
inverse of the ratio of the diffusion coefficient to the maxi- 
mal alowable, Bhom-diffusion coefficient and is always larger 
than 1. The same arguments led Berezhko & Volk (2006) 
to the conclusion that v s > 1.5 x 10 8 cm s" 1 for RX J1713.7- 
3946 (the value they obtained from the broad-band fit is 
v,w 1.8 x 10 8 cms" 1 ). 

Substituting T p = l0T p ,i keV [following Eq. (62)] in Eq. 
(12) and assuming p>2we find: 



v^Lv^ pp 



ux TB\hu=CT p 



< iorViC 1/2 rif . 



(63) 



( e is the ratio of post-shock electron and proton temperatures 
(59) and e p = 0.1e p ._i is the fraction of the thermal energy in ac- 



celerated protons. Comparing Eq. (63) to the observed ra- 
tio of fluxes per logarithmic frequency at hv = 1 TeV and 
hv = 1 keV, which for RX J1713.7-3946 and RX J0852.0- 
4622 is vf v { TeV)/V/„( keV) <~ 10 and 2 respectively, we see 
that a PP origin of the > 1 TeV is unlikely for RX J1713.7- 
3946 and not possible for RX J0852.0-4622 (since a TB flux 
greatly exceeding the observed X-ray flux would be implied) 
as long as there is significant collisionless electron heating 

If there is no collisionless electron heating, we can use Eq. 
(15) (with T p = 10r„ i keV and n = 0.1n_i cm" 3 ) 



pp 



VxL Vx TB\hv=T e 



< 0.04e ; ,,_ 1 r / , 4 f (A c/ ,, 1 . 5 «-if k yr)" 1/5 . 



(64) 



The electron temperature will be [see Eq. (14)]: 

T e > QMKxsn-iTp.if' 5 keV. (65) 

For a temperature of T e > 0.6 keV the thermal emission im- 
plied from Eq. (64) for a PP model, would likely be detectable 
in RX J 17 13.7-3946, especially if we take into account that 
there would be line emissions that would have higher lumi- 
nosities. In RX J0852. 0-4622, emission at frequencies below 



10 



1 keV might be hard to detect due to the high background 
of thermal radiation coming from the Vela SNR (Slane et al. 
2001). We note that if the proton acceleration is very efficient 
e p <~ 1, there is no collisionless heating and the TB emission 
is not considerably lower than the observed non-thermal X- 
rays, a PP origin cannot be ruled out based on this argument 
alone for either SNR. We conclude that the > 1 TeV photons 
from RX J17 13.7-3946 and RX J0852.0-4622 are unlikely to 
be emitted by PP interactions and thus are likely emitted by 
IC scattering. 

5.3. IC scenario 

We next ask whether the broad-band spectrum of these 
SNRs is consistent with an IC source of the 7-rays. 

First note that for both SNRs, the inferred cutoff in the syn- 
chrotron at <~ 1 keV is consistent with the cutoff observed in 
the <~ 10 TeV emission (we should note that for RX J0852.0- 
4622 there is only a sign of a cutoff, the uncertanties do not 
alow a firm conclusion) if we assume a magnetic field of order 
10 [iG (Aharonian et al. 2006, 2007; Porter et al. 2006, some- 
what less for RX J0852. 0-4622). In particular Eq. (47) is sat- 
isfied (up to the uncertainties in the cutoff frequencies). This 
by itself can be considered as an indication of an IC source. 

As the 7-ray observations extend somewhat below the cut- 
off, down to w 0.3 TeV, it is reasonable to compare the 
gamma ray emission directly with the radio emission, ignor- 
ing the possible suppression of the gamma ray flux due to 
cooling. Comparing equation (18) with the observed ratio 
of 7-ray to radio flux, vf v { TeN)/vf v ( GHz) ~ 100, we see 
that the expected ratio (for p = 2) is 5 - 10 times larger than 
observed in these SNRs (larger values corresponding to RX 
J0852. 0-4622). This apparent discrepancy can be due to cool- 
ing suppression of the IC flux or due to a value of p slightly 
larger than 2 (e.g. p = 2.2 would result in a factor of 5) consis- 
tent with the assumptions made here (a lower observed ratio 
would be inconsistent). 

We next note that for n ~ 0.1 cm -3 ,/? ~ 10 pc and£ ~ 10 51 
the expected cutoffs in the radio and 7-ray spectrum, Eq. (37), 
(38) are consistent with the observed cutoffs and cooling may 
or may not be important. The expected radio flux according to 
(19) is consistent with the observed ~ 1 GHz flux for the cor- 
responding distance d <~ 1 kpc. We would like to emphasize 
that there are more free parameters than constraints and these 
values are not the only ones allowable by these constraints. 

We conclude that the PP contribution to the > 1 TeV flux 
is negligible and that an IC source for the > 1 TeV flux is 
consistent with the observed broad-band spectrum. 

Using Eqs. (58) and (63), we see that the expected neutrino 
flux (being roughly equal to the 7-ray flux) is constrained for 
these SNRs to values 

£vfs„ < 10" 12 erg cm" 2 s" 1 . (66) 

The neutrino detection rate per logarithmic neutrino energy 
by a neutrino detector with an area A = A km 2 km 2 is given by 

de v 

~ 0-2 -7 — p5 : 7^km 2 Vr i 

10" 12 erg s" 1 cm" 2 km J 

(67) 

where P v ^ water is the probability that a neutrino will interact 
with the water and produce a muon within a distance from 



the detector that is smaller than the muon cooling distance, 
and is approximately given by: P VIXj water ~ 10~ 6 eJ, TeV . This 
flux is probably too low to be detected by current and planned 
neutrino observatories. 

5.4. Comparison with previous studies 

Next we compare the results presented in § 5.2 and § 5.3 to 
previous studies of these SNRs. 

In § 2.4 it was shown that Eqs. (7) and (16) agree with 
the results of studies of RX J1713.7-3946 (Aharonian et al. 
2006; Berezhko & Volk 2006; Porter et al. 2006; Moraitis 
& Mastichiadis 2006) and RX J0852.0-4622 (Enomoto et al. 
2006; Aharonian et al. 2007) to within a factor of ~ 2. We 
note that all models in which the 7-ray emission is dominated 
by PP, avoided the implied excessive synchrotron emission 
(but not the implied excessive thermal X-ray Bremsstrahlung 
emission, see § 5.2) by assuming an extremely low value of 
n~ l K ep , of n~ l K ep < 10~ 3 . Such low values of n~ l K ep are not 
plausible since a high density n 3> 0.1 is inconsistent with the 
lack of observed thermal X-ray emission and a low value of 
K ep < 10" 4 is inconsistent as shown in § 4. 

5.4.1. Claims against ICfor RX Jl 713. 7-3946 and RX 
J0852. 0-4622 

We next discuss the main claims that were raised against an 
IC source for the gamma ray emission in /?X/1713.7-3946 
and RXJ0852.0- 4622. 

Low magnetic field: As discussed in § 5.3, a magnetic field 
of B <~ 10 /xG is implied if the gamma ray emission is due to 
IC. The value of the magnetic field was estimated to be much 
higher, of order 100 /xG (Berezhko & Volk 2006; Volk et 
al. 2005; Bamba et al. 2005b) by interpreting thin filaments 
observed in the X-ray images as the result of small cooling 
lengths of the emitting electrons. If true this would rule out 
IC as the source of the gamma ray emission. The thin fila- 
ments could alternatively be interpret as thin regions of en- 
hanced magnetic field (e.g. Pohl et al. 2005) in which case 
the magnetic field cannot be estimated directly. One way to 
distinguish between the interpretations is by comparing high 
resolution radio and X-ray images. The interpretation that the 
filaments are due to cooling of multi- TeV electrons implies 
that similar features should not be seen in the radio image 
since the electrons responsible for the radio emission hardly 
suffer from radiative cooling (Vink & Laming 2003). An X- 
ray to radio comparison was done in Lazendic et al. (2004), 
however the low resolution radio images do not allow a de- 
cisive conclusion. We should note that the same arguments 
were used to deduce high magnetic fields in Tycho's SNR 
and the remnant of SN1006 (for which a high resolution ra- 
dio image exists, Dickel et al. 1991; Rothenflug et al. 2004), 
while some of the thin filaments in the X-ray emission are 
clearly seen also in the radio images , [compare (Dickel et al. 
1991) figure 1 to (Bamba et al. 2005) figure 1, and (Rothen- 
flug et al. 2004) figure 1 to (Bamba et al. 2003) figure 1, see 
also (Cassam-Chenai et al. 2007)], a fact that was ignored by 
Bamba et al. (2005) and Volk et al. (2005). 

Detailed spectral shape: Berezhko & Volk (2006) claim 
that the observed X-ray flux cannot be properly fitted for 
a magnetic field of the order <~ 10 /iG. Aharonian et al. 
(2006) claim that the shape of the gamma ray spectrum in 
RX J 17 13.7-3946 does not coincide with IC since an elec- 
tron spectrum chosen to fit the radio and X-ray observations 
produces a narrow peak in v f v in disagreement with the flat 
gamma ray spectrum observed. We do not see these claims 
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as an inconsistency as the physics of the cutoff in the particle 
spectrum is not really known. For example, the assumed dif- 
fusion coefficient value is not known for all energies. If the 
magnetic field disturbances are generated by the accelerated 
particles, the spectrum at scales relevant to the particles with 
energies close to the cutoff scale may be different than for in- 
termediate scales. In addition, if the high energy end of the 
electron energy distribution is affected by synchrotron energy 
losses, a flat vf v oc v {) IC spectrum would be expected. 

We note that if the synchrotron peak was resolved, a more 
trustable comparison of IC and synchrotron spectrums could 
have been done as long as the effect of the interstellar infra- 
red radiation is negligible (see Porter et al. 2006). 

We also note that there is some inconsistency in the model 
parameters assumed by Aharonian et al. (2006). They assume 
an age of 1000 yr, a distance of 1 kpc and an ambient density 
of n = 1 cm" 3 . For such a distance and density, the swept up 
mass is M ss (4ir/3)R 3 nm p <~ 100M Q which is clearly in the 
ST regime and implies an energy in the swept up material of 
E w 0.5nm p R 5 t~ 2 w 10 52 erg, which is rather large. Demand- 
ing an energy of 10 51 erg, for example, would imply an age of 
about 3000 yr, for which cooling in a magnetic field of 10 ^G 
may be important (the effect of cooling would be to flatten the 
IC and synchrotron peaks). 

6. DISCUSSION 

In this paper we derived simple analytic tools for analyz- 
ing the radio, X-ray and > 1 TeV 7-ray continuum emission 
mechanisms in shell-type SNRs. The emission mechanisms 
considered were Synchrotron, IC of CMB photons by accel- 
erated electrons, proton-proton collisions of accelerated pro- 
tons with ambient protons and thermal-Bremsstrahlung. In 
§ 2 we wrote down the luminosity ratios of these emission 
mechanisms (ignoring the energy cutoffs), Eqs. (7), (12), (16) 
and (18). These ratios are independent of the SNR energy 
and of the distance to the SNR. In § 3 we wrote down the 
(energy and distance dependent) expected cutoffs in the non- 
thermal radiation spectra, Eqs. (37)-(39), due to cooling and 
limited SNR age assuming DSA as the acceleration mecha- 
nism and Sedov-Taylor evolution. In addition we obtained an 
energy and distance independent constraint, Eq. (50), for the 
PP flux in case the IC spectrum is suppressed, and an energy 
and distance independent lower limit for the proton tempera- 
ture T p for SNRs in which non-thermal X-rays are observed, 
Eq.(43). We note that the synchrotron cutoff due to cool- 
ing given by Eq. (33) (a similar expression was derived by 
Berezhko & Volk 2004) naturally explains the fact that syn- 
chrotron emission does not extend to photon energies greatly 
exceeding <~ keV in known SNRs (see e.g. Reynolds & Keo- 
hane 1999; Hendrick & Reynolds 2001). This is simply be- 
cause the shock velocities in SNRs do not greatly exceed a 
few thousand km s" 1 . 

In § 4 we derived a lower limit to the value of K ep , the 
ratio of the number of accelerated electrons and the number of 
accelerated protons at a given energy, K ep > 10" 3 , by studying 
the radio observations of SNRs in M33. Here we assumed that 
the value of K ep (including the possible contributions from the 
ISM CRs) does not vary considerably between SNRs. This 
parameter enters into the ratios between IC and Synchrotron 
emissions to PP emissions. 

In section § 5 we applied the results of the earlier sections 
to show that the broad-band spectra of the SNRs, RX J 17 13.7- 
3946 and RX J0852.0-4622 are inconsistent with a PP origin 
and are consistent with an IC origin of the > TeV emission. A 



PP dominated TeV emission would imply radio synchrotron 
and probably thermal X-ray Bremsstrahlung fluxes that would 
greatly exceed the observed X-ray flux. 

The neutrino flux from these SNRs is expected to be lower 
than e v f £v < 10" 12 erg cm" 2 s" 1 and is probably too low to be 
detected by current and planned neutrino observatories. 

We compared our main results with previous studies of 
these SNRs (Tables 1 and 2) and showed that our simple 
analytical expressions are in good agreement with more de- 
tailed calculations. All models, in which the 7-ray emis- 
sion is dominated by PP, avoided the implied excessive syn- 
chrotron emission (but not the implied excessive thermal X- 
ray Bremsstrahlung emission, see § 5.2) by assuming an ex- 
tremely low value of n~ l K ep , n~ l K ep < 10" . Such low values 
of n~ l K ep are not plausible since a high density n ^> 0.1 is in- 
consistent with the lack of observed thermal X-ray emission 
and a low value of K ep < 10~ 4 is inconsistent with radio obser- 
vations of SNRs in nearby galaxies as shown in § 4. Previous 
claims, that the 7-ray emission in SNRs RX J 17 13.7-3946 and 
RX J0852.0-4622 is not IC where discussed in § 5.4.1. 

Interpretation of the narrow filaments seen in the X-ray pic- 
tures as cooling width of the emitting electrons was used to 
obtain magnetic field estimates of order ^100 [iG (Berezhko 
& Volk 2006; Volk et al. 2005; Bamba et al. 2005b). This 
would rule out an IC source and thus seem implausible. As an 
illustration, this would require a value of the electron:proton 

ratio of K ep <~ 10" 5 «_iB_^ 2 to explain the <~ 100 ratio of TeV 
to GHz fluxes per logarithmic frequency (without solving the 
thermal-Bremsstrahlung problem). The interpretation of the 
narrow filaments as cooling width of the multi-TeV X-ray 
emitting electrons implies that similar filaments are not ex- 
pected in the radio observations. There are at least two exam- 
ples (Tycho's SNR and the remnant of SN1006) where similar 
filaments are observed in both radio and X-rays. This puts into 
question the high B interpretation of the X-ray filaments (see 
discussion in § 5.4.1). 

We note that the Synchrotron to PP ratios would be affected 
if the magnetic field is enhanced in a small region behind the 
shock as suggested above but that the conclusion that a PP 
model requires low values of K ep would not change. To see 
the effect of thin enhancement regions, assume an extreme 
case where there is a strong magnetic field B in a small re- 
gion d <C behind the shock, and a negligible magnetic field 
elsewhere. Assuming that the accelerated electrons are not 
confined to this region, the radio emission would be propor- 
tional to dB 3 > 2 and it would be possible to allow for a higher 
value of K ep in a PP model for a given value of the magnetic 
field. Note however, that in order to cool the electrons emit- 
ting the TeV IC for a given SNR age (see discussion in § 3.4), 
the magnetic field would have to be larger in order to cool the 
electrons in the time they reside in the high magnetic field re- 
gion and thus will have to be larger by a factor oc c/" 1 / 2 com- 
pared to a homogenous case. So the suppression of the radio 
flux due to the small emitting region, given that the IC emit- 
ting electrons are cooled, will be roughly equal to (d/R)~ x l A 
where R is the remnant radius. The thin filaments observed 
by Chandra, have widths of 2' and 1' for RX J 17 13.7-3946 
and RX J0852.0-4622 respectively (Berezhko & Volk 2006; 
Volk et al. 2005; Bamba et al. 2005b). Taking into account a 
projection factor of w 7 (Berezhko & Volk 2006) the emis- 
sion regions widths are fractions d/R <~ 10" 2 and 3 x 10" 3 of 
the SNR radii respectively. This would require a correction 
factor of {d/R) l l A ^3-5 to equation Eq. (59) and will not 
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change the conclusions. Furthermore, assuming that the X- 
ray synchrotron emitting electrons are also effectively cooled 
in this region, Eq. (60) will remain valid. 

Using the magnetic field value B ~ 10 yuG, the ratio of mag- 
netic field energy to thermal energy of swept up material is 
roughly given by: 



e B „ 3^-R 3 E^ ept ~4x 10~ 4 B 2 _ 5 R]e;1 



swept.5 1 1 



(68) 



where Ewept = 10 51 £ swep t,5i erg is the total energy in swept 
up material. We note that if the distances to these SNRs 
are a few kpcs e B would equal a few percents. For exam- 
ple, a radius of R = 2>QR\ 5 pc, implying distances of 3R\ 5 kpc 
and 1.5/?i. 5 kpc to RX J17 13.7-3946 and RX J0852.0-4622 
respectively, implies e B <~ 0.01B^ 5 /?j 5 £^ ept51 and is consis- 
tent with all observations. We note that larger distances imply 
smaller densities since the velocity is limited from below by 

Eq. (61), Vj > 3 x 10 8 i^y cm s" 1 , and the number density can 
roughly be expressed as n <~ 2 x 10~ 3 £ swept 5i Vg^m" 1 /?"^ cm" 3 
where v, = 3000vg 5 km s -1 . Such low densities are expected 
if these shocks are propagating into progenitor winds (see e.g. 
Berezhko & Volk 2006, and refferences within). 

7-ray observations in the GeV to sub-TeV range by the 
GLAST experiment will hopefully allow a clear direct distinc- 
tion between the IC predicted spectrum, vf v (XV 1 / 2 (which 
is thus predicted for the SNRs RX J 17 13.7-3946 and RX 
J0852.0-4622) and the PP predicted spectrum vf v oc v° (with 
a cutoff at ~ 100 MeV energies). Using Eq. (7), the expected 
IC to PP flux ratio for GeV photon energies is approximately 



L v i C /L v PP ( GeV) w 3K ep - 2 n_ x . 



(69) 

A non-negligible contribution of the PP emission cannot be 
ruled out (for smaller photon energies the PP emission is 
strongly suppressed). However, it is certainly possible that 
PP emission is masked out by IC at all photon energies for 
these SNRs. 

An interesting question is what kind of SNR parameters 
are required in order to have an observable gamma ray emis- 
sion dominated by PP collisions. Higher densities would re- 
sult in higher PP emission albeit with lower maximal pro- 
ton energy. The maximal proton energy is proportional to 
£p,max oc BRv s oc EBR~ X l 2 n~ x I 2 . Based on the observation that 



electrons are accelerated to <~ 60 TeV energies in these SNRs 
we assume protons are accelerated to similar energies (prob- 
ably somewhat higher if the electrons are limited by cool- 
ing). Therefore, comparing to these SNRs, we have freedom 
to increase the density by a factor of <~ 100 (fixing the en- 
ergy, radius and magnetic field), while keeping protons ener- 
getic enough to produce ~ 1 TeV photons. The cutoff pho- 
ton energies in the IC spectrum and the synchrotron spec- 
trum are both proportional to oc v^ZT 1 oc ER~ 3 B~ l ri~ 1 and 
oc v 2 B 2 R 2 oc ER~ l B 2 n~ l , for cooling and age limits respec- 
tively [see Eqs. (37) and (38)]. A factor of ~ 100 in the den- 
sity (for fixed energy, radius and magnetic field) would shift 
the IC and synchrotron cutoff energies by a factor of 1/100, 
strongly suppressing the TeV IC and keV synchrotron emis- 
sions. At the same time, a larger density would increase the 
thermal X-ray emission (as long as the post shock tempera- 
ture does not fall below the X-ray observable energies). It 
is therefore likely that SNRs with considerably higher ambi- 
ent densities have observable PP dominated TeV emission. 
Such SNRs will have thermal or no observable X-ray radi- 
ation rather than non-thermal X-ray radiation. At <~ 1 GeV 
photon energies, densities exceeding n > 03K~p_ 2 [cf. (69)] 
are enough for PP emission to dominate the IC emission. 

Neutrino emission from PP collisions is similarly expected 
to be higher in SNRs evolving in high density environments 
(the neutrino flux roughly equals the PP gamma ray flux) and 
are likely to be better observed in SNRs with strong thermal 
X-ray emission (or no X-ray emission). For SNRs with ob- 
served thermal X-ray emission, the expected neutrino flux can 
be estimated directly using (12). 

We conclude that there is need for a detailed analysis us- 
ing the X-ray and radio data of SNRs in order to find suitable 
candidates for PP 7-ray and neutrino emission. The analyti- 
cal tools developed in this paper may be used to estimate the 
expected 7-ray and neutrino fluxes and to determine the domi- 
nant 7-ray emission process based on existing radio and X-ray 
observations of SNRs. 
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APPENDIX 
EMISSION MECHANISMS 

Thermal Bremsstrahlung 

The thermal Bremsstrahlung emissivity per unit frequency of an optically thin plasma with temperature T e is given by (Rybicki 



& Lightman 1979): 



3m e c 3 V 3m, 



1/2 



(Al) 



where n e , T e are the electron number density and temperature respectively, n,, Z are the ions' number density and charge respec- 
tively, and gff is the thermal Gaunt factor. For a plasma consisting of electrons and protons with equal number density n we 
have 



ve f J = 



, oja.eC . 
3tt V T, 



2\ 1/2 , 

m e c z -\ hv m _ h „ lT 
n — T,ne hv ' T ' 



8ff- 



(A2) 



The function xe x attains its maximal value (e ') at x = 1 . The maximal Luminosity per logarithmic frequency is thus 

vL,f J,hu=T, = V ^e~ l a e g ff N(j T nc^ m e c 2 T e . (A3) 
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For hv = T e , 100 eV < T e < 10 keV, the value of g ff is in the range, 0.8 < g ff < 1.2 (e.g. Karzas & Latter 1961). 

We next consider the expected value of the electron temperature due to coulomb heating by protons. The equation for the 
change in the electron temperature due to coulomb collisions with protons is given by (eg. Ichimaru 2004): 

^ = ( r,-r«)^(ii + i)" M A„, ,A4) 

at 3m e m p \m e m p J 

where X ep is the Coulomb logarithm. Assuming that m e /m p <C T e /T p <C 1 the electron temperature after a time t = ?k yr kyr will 
be: 

T e ~ 0.6 (A Vi i.5nofkyr7>.kev) keV, (A5) 

where X ep = 30\ ePt \.5. For 10" 2 < T e /T p < 0.6 the correction to this expression is smaller than 20%. 
Next consider power law distributions of accelerated electrons or protons, 

dNi A -P i A 

-I- = A Hi , (A6) 

where i = e,p. 



Gamma rays from proton-proton collisions 
The spectrum of emitted photons is given by Drury et al. (1994) : 

1 i I \P inel dN p . 

-L„ P p = (mx)' 7 a pp nc -^r \hu, 



K^^h;, (A7) 



where a l pp is the inelastic proton-proton cross section and {mx) P s is the spectrum weighted moment for particles of type S. This 
can be written as: 

vh v p P = C PP (p)2s p ^a i p fnch^, (A8) 

CIS p 

where e p dN p /de p is to be evaluated at e p (v) = lOhv. Ignoring the correction factor Cpp(p), this is equivalent to assuming that 
each inelastic p-p collision produces two photons with energy hv = e p /\0. The correction factor is given by 

For p = 2,2.2 we have C PP (2) w 0.85,C PP (2.2) w 0.66 (values of (mx)^ were taken from Drury et al. 1994). 

IC radiaion ofCMB photons 

The spectrum of IC scattered photons of a black body target with temperature T in the Thompson regime is given by (Rybicki 
&Lightman 1979) 

\h v ,c=A e ^(Tf%(p)(^ 1)/2 , (A10) 

where 

f(p) = 2^ p2 : 4p+U r (Z±*) ( ( p -^) . (Al 1) 

KU (p + 3)\p + 5)(p+l) \ 2 2 J 

1 dN e 4 , 

vU ic = Ci C {p)-e e {v) — -a Tl 2 e {v)U T c, (A12) 

2 de e 3 

where e e dN e /de e is to be evaluated at e e {v) = 7 e (i/)m e c 2 = m e c 2 (hv /3T) X I 2 and Ut = aT A is the energy density in the black body 
photons. Ignoring the correction factor Cic(p), this is equivalent to assuming that each electron emits all the power ^ctt^Utc, 
in photons of energy hv = j 2 3T. The correction factor is given by: 

C 1 c(p) = 3- {p - 1) / 2 ^ir- 4 F(p). (A13) 

For 2 < p < 2.2 we have Qc(p) ~ 0.8 (to within 5%). It is useful to note that ^ 2 {v)Uj = [Ut /(3T)]v w 0.9n T hv where n T is the 
number density of black body photons. 



This can be written as, 
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Synchrotron radiation 
The spectrum of synchrotron radiation is (Rybicki & Lightman 1979) 



1 



. p 19 \ /p 1 ^ \/3? 3 BsinQ: f lixmcv 



p+\ V 4 12 



T2) 



3qB sin a 



-(p-D/2 



where a is the angle between the electrons velocity and the magnetic field direction. Using, 



1 f\ . ... -ci V^T(^) 
— / 27rsinaaa:(sina) 2 = — , 



the spectrum for an isotropic distribution is 



L v =A e D{p) 



q 3 B ( mcv\ 
mc 2 \ qB J 



-(p-l)/2 



with 



D(p) = 



1 V^(^)r(f+B)r(f-£)^ 27r 



p+i 



2 r(^) 



-(p-D/2 



This can be written as 



1 dN e 4 , 

vL v = C S yn(p)~e e (is)—-aT%(is)U B c, 

2 de e 3 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



where e e dN e /e e is to be evaluated at e e (v) = "/ e (v)m e c 2 = (2v/v B ) l / 2 m e c 2 , v B = qB/(2TTtn e c) and Ub = B 2 /(8tt). Ignoring the 
correction factor Cs yn (p), this is equivalent to assuming that each electron emits all it's power of (4/3)arj 2 U B c in photons of 
energy hv = r / 2 hv B j2. The correction factor is given by 



C Sya ( P )=^x(4w)^D(p) 



(A19) 



and is approximately Cs yn (/?) ~ 0.8 (to within 5%) for 2 < p < 2.2. 
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